METHOD AND SYSTEM OF DEFINING A MODEL OF ONE OR MORE 

ORGANS 



FIELD OF INVENTION 

5 

The invention relates to a method and system involving defining a model of one or more 
organs or part(s) thereof then using the model to interpret remote electrical and/or 
magnetic recordings. The invention has primary application in data interpretation 
relating to the gastrointestinal tract. The invention is equally applicable to any other 
10 organ from which magnetic or electrical activity can be recorded, for example the heart, 
brain and uterus. 

BACKGROUND TO INVENTION 

15 Prior art techniques exist to obtain electrical and/or magnetic recordings from gastric 
activity, from intestine activity, from cardiac activity, and many other musculature 
organs. 

The twelve lead ECG (electrocardiogram) is used to measure cardiac electrical activity 
20 and is a standard that has been almost universally adopted. At present, there is no 
standard for measuring an EGG (electrogastrogram). The low conductivity tissues that 
lie between the active tissue in the gut and any cutaneous electrodes hinders the recording 
and subsequent interpretation of EGG signals. The magnetogastrogram (MGG) does not 
suffer from these problems, as the permeability of biological tissues is very similar to that 
25 of free space. 

It would be particularly advantageous to create a computer model that is capable of 
reproducing and/or interpreting MGGs and electrical and/or magnetic recordings of 
various other body organs. 

30 



2 



SUMMARY OF INVENTION 

In broad terms in one form the invention comprises a method of defining a model of one 
or more organs or part(s) thereof from multiple images of the organ(s) or part(s) thereof, 
5 the method comprising the steps of generating a computational mesh of one or more 
organs from multiple images of the organ(s); generating a representation of musculature 
associated with the organ(s); calculating electric and/or magnetic fields associated with 
the muscle layers; and defining a model based on the computational mesh, and the 
electric and/or magnetic fields. 

10 

The invention also comprises a method of estimating the location of one or more sources 
of magnetic and/or electric fields in a subject by defining a model of one or more organs 
or part(s) thereof as described above, obtaining one or more measured magnetic and/or 
electric fields from a subject, and estimating the location of one or more sources of 
15 magnetic and/or electric fields based at least partly on the model of one or more organs 
and the measured magnetic and/or electric fields. 

In broad terms in another form the invention comprises a model defining system for 
defining a model of one or more organs or part(s) thereof from multiple images of the 

20 organ(s) or part(s) thereof, the system comprising a mesh generation component 
configured to generate a computational mesh of one or more organs from multiple images 
of the organ(s) and a representation of muscle layers associated with the organ(s); an 
electric/magnetic field component configured to calculate electric and/or magnetic fields 
associated with the muscle layers; and a model creation component configured to define a 

25 model based at least partly on the computational mesh and the electric and/or magnetic 
fields. 

In another form the invention comprises a source location system for estimating the 
location of one or more sources of magnetic and/or electric fields in a subject. The 
30 system comprises a model defining system as described above and a location estimator 
configured to estimate the location of one or more sources of magnetic and/or electric 
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fields based at least partly on the model of one or more organs and data obtained from 
one or more measured sources of magnetic and/or electric fields from a subject. 

BRIEF DESCRIPTION OF THE FIGURES 

5 

Preferred forms of the invention will now be described, by way of example, with 
reference to the accompanying figures in which: 

Figure 1 shows a preferred form method for defining a model in accordance with the 
10 invention; 

Figure 2 shows a preferred form system in accordance with the invention; 
Figure 3 illustrates an image of one or more organs of a subject; 

15 

Figure 4 shows a preferred form fitted stomach mesh; 
Figure 5 shows the digitised and fitted outer skin surface of the torso; 
20 Figure 6 illustrates individual meshes of multiple organs of the digestive system; 
Figure 7 illustrates a representation of musculature; 
Figure 8 shows an example of a high resolution finite difference mesh; 

25 

Figure 9 illustrates the results of a simulation in accordance with the invention; 
Figure 10 illustrates magnetic field vectors generated in accordance with the invention; 
30 Figure 1 1 illustrates point source localisation in accordance with the invention; 
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Figure 12 shows a sample output generated in accordance with the invention; and 

Figure 13 shows a further sample output generated in accordance with the invention. 
5 DETAILED DESCRIPTION OF PREFERRED FORMS 

Figure 1 sets out a preferred form method of the invention for defining a model of one or 
more organs or parts of organs of a subject. The first step is optionally to complete a 
subject record 10 containing as much relevant information as possible, including weight, 
10 height, age and current medication. The subject record could also include any history of 
problems associated with the organ(s) of interest, for example digestive problems. 

The next step is to obtain images of the subject 20 to obtain the anatomical and 
geometrical information that is necessary to create a subject-specific model of the 
15 digestive or other system and surrounding tissues. The imaging modality used to gain 
this information on the internal composition of a subject includes one or more of MRI 
(magnetic resonance imaging), CT (computer tomography), ultrasound, PET (positron 
emission tomography) and X-ray imaging. 

20 It is also envisaged that the data obtained from the selected imaging modality could be 
supplemented in one or more ways. For example, an additional imaging technique could 
be used such as a laser scanner or camera system to capture the external surface of the 
subject. Distance measurements could be manually taken from the subject, such as 
lengths and circumferences. A mechanical or other acquisition device could be used to 

25 record the location of fiducial markers to aid in the registration of gathered data. 
Furthermore, the body fat of a subject could be measured. 

A contrast agent could be used when imaging the subject such as a barium swallow to 
better elucidate the geometry of the digestive system, or other technique suitable for the 
30 organ(s) of interest. 
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Following the imaging process, the outputs from the imaging of the subject are likely to 
be either 2-dimensional image slices that are not restricted to the same orientation, or 3- 
dimensional volume image sets. These outputs could display one or more organs or one 
or more parts of organs of the subject imaged. 

5 

The information within these images that is of interest for one application of this 
invention is the location and anatomical configuration of the organs that make up the 
digestive system. Also of interest are the organs that surround the digestive system and 
influence the shape or functional properties of the digestive organs. 

10 

For each organ or part of interest, representations of the surfaces of the organ are 
required, as well as any internal surfaces that are required to represent regions with 
different physiological properties. It is preferable to represent at least the electrically 
active organs of the digestive system of a subject and all of the surrounding organs that 
15 influence magnetic fields and electric fields associated with these electrically active 
organs. The electrically active digestive organs include the oesophagus or esophagus, the 
stomach, the small intestine (duodenum, jejunum, ileum) and the large intestine (colon). 
The tissues of interest surrounding these organs include the liver, the pancreas, the 
abdominal muscles, the subcutaneous fat, and the fat within the abdominal cavity. 

20 

The next step is to generate 30 a computational mesh representing the above features. 
The generation of the computational mesh is further described below. There are several 
different ways of creating the computational mesh from multiple images, for example 
image slices, of the organ or organs. 

25 

One such method is to automatically generate computational meshes using software 
which automatically segments image data and an automated mesh creation process that is 
subject to constraints depending on the precision requirements of the mesh. 

30 Another method is to create data points from an image set and then fit an initial generic 
model to the data, minimising the distance between the data and the fitted model. These 



65295-1 



6 

data points can be created manually using digitising software, for example the software 
described in patent specification WO 01/01859 to Auckland UniServices Limited. It is 
envisaged that if the gathered data is too sparse in some areas, perhaps because of an 
imaging artifact, it is possible to perform a second type of fit known as anthropomorphic 
5 fitting. This type of fitting uses a set of fiducial markers to adjust a previously fitted 
organ mesh so that it matches the current geometry. 

In each case, a mesh generation component could be a computer program installed and 
operating in computer memory which is configured to generate a computational mesh of 
10 one or more organs from multiple images of the organ(s). 

Following generation of the computational mesh for organs, or parts of organs, the next 
step is to generate 40 a representation of musculature as will be described below. The 
musculature generating the magnetic fields are either inside or form part of the organs. 
15 There are also muscles outside the organs that act in a passive manner, such as the 
abdominal wall muscles, but these are not a primary field source. The mesh generation 
component could also be configured to generate representations of the muscle layers 
associated with the organ(s), either inside or forming part of the organ(s). 

20 It is often difficult to image a subject in sufficient detail to obtain all the required 
information, and so it is preferable to add prior knowledge to the anatomically accurate 
geometrical model. This information could include data on the cellular and tissue 
composition of each organ and the spatial and temporal variations that occur. In 
muscular regions, the microstructure of the muscles needs to be added as preferential 

25 directions of electrical propagation and contraction are often present. 

A combination of accurate geometry and the appropriate physiology contributes to the 
generation of a specific subject model of the system of interest. This in turn provides a 
framework in which the equations governing the processes representing gastric activity, 
30 intestinal activity, cardiac activity, or other activity can be solved accurately. 
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Referring to Figure 1, the next step is to calculate 50 the magnetic fields associated with 
the model layers generated in steps 30 and 40 and to define 60 a model based on the 
computational mesh, the representation of musculature, and the electric and magnetic 
fields. 

5 

Figure 2 illustrates a preferred form system 200 for carrying out the method described 
above. Images 202 represent, for example, multiple image slices of an organ or organs. 
Images 202 and optionally prior knowledge data 204 is input to mesh generation 
component 206. The mesh generation component 206 is configured to generate a 
10 computational mesh 210 and optionally representations of the muscle layers 212. 

As will be described below, it is envisaged that data obtained from non-invasive 
electric/magnetic measurement data 214 be used as input to an electric/magnetic field 
component 216. The electric/magnetic field component 216 could include a computer 
15 program in which mathematical models and related mathematical equations are 
calculated. The component 216 is configured to calculate electric and/or magnetic fields 
218 associated with muscle layers. 

The computational mesh 210, the electric and/or magnetic fields 218 and optionally the 
20 muscle layer representations 212 are then input to a model creation component 220. The 
model creation component is preferably a software-based component configured to define 
a model based on the above inputs. Model 222 preferably represents one or more organs 
or part or parts thereof. 

25 The system 200 optionally further includes a location estimator 224. The estimator 224 
preferably comprises a computer program that is configured to solve an inverse problem. 
The estimator 224 takes as input the model 222 generated by component 220 and 
electric/magnetic measurement data 214 and estimates the location of one or more 
sources of magnetic and/or electric fields shown in Figure 2 as electric/magnetic field 

30 location(s) 226. 
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Figure 3 illustrates an image 100 of one or more organs of a subject. Photographic 
images in the axial plane of a subject that are available at 1mm intervals over the length 
of the subject body have been digitised. The components of interest are traced on every 
second transverse digital human image which translates to a vertical resolution of 2mm as 
5 the original slices were taken at 1mm intervals. 

The components particularly of interest in gastrointestinal activity are the oesophagus, the 
stomach, the duodenum, the jejunum and the ileum. The outer surface of the oesophagus, 
stomach, duodenum and small intestine are traced, and the centre line of the remainder of 

10 the small intestine is also located. It will be appreciated that the components of interest 

1 will vary depending on the activity under analysis. 

For the oesophagus, stomach and duodenum, initial linear quadrilateral surface elements 
are created by selecting a regular array of data points to be the initial nodal positions. 
15 From this, the nodal positions and then nodal derivatives are fitted to the digitised data set 
for each component to create a bi-cubic Hermite outer surface description of the system. 

The RMS (root mean square) errors between the final fitted surfaces and the digitised 
data are less than 1mm for each of the meshes. From the outer surface, a volume mesh is 
20 created through an inward cylindrical projection and a wall thickness of 5mm is chosen 
as an average value. 

In Figure 3 the outer wall of the stomach has been digitised. Graphics window display 
300 shows the currently digitised points in 3-dimensions and provides an indication of 
25 where data may be sparse, missing or incorrectly placed. An image control dialog box 
310 enables the manipulation of an image set and also allows multiple organs to be 
digitised into separate data groups on a single image. Once the data has been collected, 
the initial computational mesh is generated to which the data is fitted. 

30 Figure 4 illustrates at 400 a preferred form fitted stomach computational mesh which is 
the result of the fitting procedure to the outer wall of the stomach. A C 1 continuous 
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quadrilateral description of the surface is calculated which is then used to generate each 
of the individual points visible on the surface. The average error between the data points 
and the surface projection of these data points on the mesh is approximately 0.6mm, 1 

5 A volume computational mesh can also be generated by either applying the same 
digitising and fitting process to the inner surface of the stomach, or by projecting the 
outer surface inwards based on known information about the thickness of the stomach 
wall. 

10 It is also envisaged that the outer skin surface of the torso be digitised and fitted as shown 
in Figure 5. As all the data is referenced to a common origin, the position of the organs 
within the torso is consistent with the original images. 

As shown in Figure 6, individual meshes of multiple organs of the digestive system can 
15 be generated and combined with a mesh of the outer skin surface of the torso. The 
combined computational mesh 600 could include the oesophagus 602, the stomach 604, 
the duodenum, jejunum and ileum making up the small intestine 606 and the colon or 
large intestine 608. The sigmoid colon has been omitted from this image for clarity. 

20 The centre line of the jejunum and ileum is used to create a topologically cylindrical 
mesh through a radial projection. The outer diameter of the jejunum is 40mm with a wall 
thickness of 4mm and the diameter of the ileum is 37.5mm with a wall thickness of 3mm. 
The transition from the jejunum to the ileum is performed on a length basis as the 
jejunum occupies approximately 40% of the length of the small intestine whereas the 

25 ileum occupies approximately 60% of the length. 

It is envisaged that all the above mesh components can be customised through a host 
mesh fitting technique to individual subjects to account for subject variability. 
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Figure 7 shows a representation 700 or model of the musculature in a small section of the 
stomach wall. The representation shows longitudinal muscle layer 702, circular muscle 
layers 704 and 706, and interstitial cells of Cajal 708 and 710. 

5 Once the combined computational mesh is defined, and representation of musculature 
generated, it is then necessary to determine the electrical activity within the organ(s) or 
part(s) thereof, for example the stomach and intestine and the related muscle layers, 
which is producing a known or measured magnetic or potential field on or near the torso 
surface. Sensors can be placed near or on the torso to obtain non-invasive measurements 
10 about the electrical or magnetic fields generated from within the torso. Electrical fields 
can be measured in a non-invasive manner by electrodes placed directly on the skin 
surface. 

Known electrical mapping systems are able to record from one to many hundred 
15 electrodes simultaneously. These mapping systems have been typically designed for 
recording electrical activity from the heart or brain. These measurements are obtained 
relative to a single or combination of electrodes recorded at the same time. 

Magnetic sensors are capable of recording changes in magnetic fields without direct 
20 contact with the torso surface. These fields can be recorded using recording devices 
known as SQUIDS, for example as described in US patent specification 5,771,894 to 
Richards et al. 

The locations on which these electrodes or sensors are positioned relative to the torso are 
25 also critical for subsequent measurement. As electrodes are placed directly on the skin 
surface, the relative positions of these electrodes will shift for each subject. On the other 
hand magnetic sensors are non-contact and are typically in a fixed position relative to 
each other. This means that their relative positions only need to be determined once but 
their overall positions relative to the torso must be determined each time. 

30 
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The EGG (electrogastrogram) and the MGG (magnetogastrogram) record not only the 
surface projections of the summated action potentials of gastric smooth muscles, but also 
a compound signal from different electrical sources in the torso. These magnetic and 
potential fields are attenuated and filtered by the surrounding organs, for example muscle 
5 and fat layers. The field data can be interpreted by pattern matching against known 
signal sets, or by using a mathematical model to interpret the data directly. 

The process of pattern matching, or comparing signal traces between a health normal and 
an abnormal patient is commonly used to determine the presence of an abnormality. For 

10 such a process to be a success, the person interpreting the data needs to take into account 
many factors, including the size of the subject, location and relative placement of sensors 
or electrodes, age and so on. It is assumed that the sensors or electrodes are placed at 
standard locations so that results can be compared between patients. The process can 
often be flawed and require significant training and an established and consistent 

15 database for comparison of results. 

An alternative method of interpreting the data is through use of a mathematical model 
where measured data can be used to directly interpret the field data. The use of a 
mathematical model makes such an interpretation less subjective to human judgement, as 

20 results are typically displayed in a more direct and meaningful manner closely related to 
the underlying origins of the events. As described above, electric/magnetic field 
component could include a computer program in which the mathematical model and 
related mathematical equations are calculated. The model creation component is 
configured to define a model based on the computational mesh, the magnetic fields 

25 calculated by the following equations and/or the representation of muscle layers. 

Assumptions are made as to the equations that govern the electric and magnetic activity 
in the subject and these equations are then solved once the computational mesh has been 
created. 

30 
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The equations that govern all magnetic and electric fields are known as Maxwell's 
equations. The frequencies of the biological signals that are of interest are generally less 
than 100 Hz, and the magnetic permeability of biological tissue is very similar to that of 
air. On this basis, Maxwell's equations may be simplified to what is termed the quasi- 
5 static assumption. 

The governing equations in differential form can be written as: 

V • E = — (1) 
10 V x E = 0 (2) 

V • B = 0 (3) 
VxB=// 0 J (4) 

15 

E is the electric field intensity, B is the magnetic flux density, J is the electric current 
density, p is the electric charge density, eo is the permittivity of free space, and (io is the 
permeability of free space. 

20 In addition to the above equations, the continuity equation is needed to ensure there is no 
build up of electric charge within a region: 



V-J = 0 (5) 

25 Three main formulations arise from the above equations. The first formulation is a set of 
equations known as the bidomain equations that govern the spread of electrical activity 
within excitable tissue. The second formulation is a generalised Laplace equation that 
describes the current flows within passive tissue regions. The third formulation is an 
equation which takes the electric fields as input and then calculates the magnetic field 

30 generated by the electrical activity. 
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The first formulation of equations model active tissue as two inter-penetrating domains 
that occupy the same physical space. The relationship between potentials in the two 
spaces across the cell membrane is the first equation in the bidomain system. 

V m =<P,-<P e (6) 

The intracellular domain represented as i represents the interior of the cells, and the 
extracellular domain represented as e represents the material surrounding the cells. Other 
terms contain m to indicate they are properties of the cell membrane. 

The two bidomain equations are: 



V-(or,VVj + V-(ff,V*.)= A fc m ^+ /,, 



ion 



+ /„ (V) 



dt 

15 V • ((a, : + a e )V <p e ) = -V (<r,.V V m )- I s2 (8) 

Here the a terms are tissue conductivities which in general will be tensors, the <p terms are 
potentials, the V m term is the transmembrane potential, the potential difference across the 
cell membrane, A m is the surface to volume ratio of the membrane and C m is the 
20 membrane capacitance. Individual cellular models are able to plug directly into these 
equations through the I lon term in the first equation. 

At a fine scale, each cellular model is able to incorporate complex subcellular processes. 
Externally applied currents may be injected into either domain through I s i or I s2 . These 
25 equations are solved using either the finite element-based finite difference method or the 
structured finite element method. 



Figure 8 shows an example of a high resolution finite difference mesh created over the 
finite elements of the stomach using the first formulation. 
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The second formulation that arises from Maxwell's equations is a generalised Laplace 
equation: 

V-( ( r„V<p„)=0 (9) 

The o subscripts denote quantities outside the active region. Bidomain equation (8) that 
solves for q> e is directly coupled to the external passive regions through the interface 
conditions on shared boundaries: 

<P* = <P 0 (10) 
a e Vp e -n e = <r 0 Vp 0 -n 0 (11) 



15 Similar formulae apply to interfaces between passive regions, ensuring that the potential 
fields and current flows are properly conserved. For electrically isotropic regions, the 
conductivity reduces to a single value and the boundary element method is used to solve 
these equations. For electrically anisotropic regions the conductivity is a 3x3 tensor and 
the finite element method is used. 



10 



20 



25 



The third formulation defines the magnetic field generated by the electrical activity and is 
defined in terms of the curl of a vector field A: 

B = V x A (12) 
The vector potential field A is then defined in terms of the electric current density: 

V 2 A = -fij (13) 



30 The current density comprises two components. The first component is the contribution 
of the primary current sources, in this case a function of the transmembrane potential 
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gradient. The second component is the contribution from distribution of the resistive 
network within the torso volume, or other relevant body volume: 

J = J p -aVp (14) 

5 

The above equations are solved using the finite element and boundary element meshes 
that are used to model the passive electric fields. 

A typical solution consists of four sets of calculations that must be performed at each 
10 time step. In the first of these calculations, the cellular terms are updated throughout the 
active region. In the second set of calculations, the transmembrane potential is calculated 
from the cellular terms and known diffusive properties. In the third set of calculations, 
the coupled extracellular/passive torso problem is solved to create a continuous electric 
field throughout the torso. In the fourth set of calculations, this electric field is used as an 
15 input to the calculation of the magnetic field. 

A simplification to this coupled system involves the introduction of equivalent dipolar 
source terms to represent the contribution of the active region. In this case, the active 
region is solved in isolation and equivalent sources are calculated through the vector 
20 summation of the primary cellular sources. Typically, this process produces in the order 
of tens of dipole sources to represent the electrical activity. These sources are then 
placed into the passive torso model and the appropriate equations are solved to obtain the 
electric and magnetic fields within and surrounding the torso. 

25 Figure 9 illustrates the results of a simulation in accordance with the invention. The bulk 
of the electrical activity is located halfway along the duodenum at the beginning of the 
small intestine within the torso. From this cellular activity, equivalent dipole sources are 
used as inputs to the passive torso solution that generates the electric fields throughout 
the torso. The primary sources and these electric fields are then used to generate the 

30 magnetic field just outside the surface of the torso. In the figure, the magnetic field is 
drawn as arrows that have been seeded at random points within the square in space. Each 
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arrow points in the direction of the calculated magnetic field and the length of the arrow 
indicates the strength of the field. 

The invention is particularly suited to solving an inverse problem. An inverse problem is 
5 a general term for a class of problems which attempt to determine the sources or events 
which produce a known result. Mathematically they can be written in the general form 
A x = b where x is an unknown variable, b is a known or measured variable and A is a 
function which is capable of mapping between the two variables. 

10 The inverse problem for the stomach and intestinal system involves determining the 
electrical activity within the stomach and intestine which is producing a known or 
measured magnetic or potential field on or near the torso surface. In the case of the 
above equation, x represents the unknown electrical sources in the stomach and/or 
intestine, b represents a known or measured resultant field (external electrical potential or 

15 magnetic fields) while A is a function which relates two fields which takes into account 
the geometries of the organs in the torso and their electrical conductivities which effect 
the path and pattern of the electrical/magnetic activity as it disperses from the site of 
activation. 

20 Inverse problems also fall into a class of problems which are known as ill-posed. This 
means that small amounts of error in the solution process and input data can result in the 
large and disproportionate errors in the computed solution. This means they are 
mathematically difficult to solve. 

25 The use of a mathematical model makes interpretation of data less subjective to human 
judgement as the results are typically displayed in a more direct and meaningful manner 
closely related to the underlying origins of the events. 

Inverse algorithms are used to compute a "source" given a known geometry and 
30 surrounding potential or magnetic field. Such a source is usually a simplified, at a spatial 
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scale larger than that of a cell, but realistic representation of the underlying events 
actually occurring in the body. 

Point source inverses locate dipoles at sites of interest. Dipole sources are a 
5 mathematical way of representing the strength and direction of a field with only a few 
parameters. A dipole is essentially a vector quantity and defined by six components. 
Three components define the centre for example (x, y, z) and three components define the 
orientation for example (dz, dy, dz) assuming a rectangular Cartesian space. 

10 Dipoles are commonly used in detecting events and localising regions of interest in the 
brain from both electrical and magnetic fields and used to a lesser extent with electrical 
fields for imaging the heart. 

A non-linear optimiser is used which minimises the difference between a 
15 known/measured potential/magnetic field and that computed from an estimated source 
applied within the torso by adjusting the source parameters. Thus, from an arbitrary 
initial estimate of a dipole(s), the resultant potential/magnetic fields are computed, and 
then sequentially adjusted until the potential/magnetic fields match. These can be 
described mathematically as: 

20 

Minimise F = f(B,B') + lambda f(Phi, Phi') (15) 
w.r.t. dipole parameters 

Where B and B' are the known and computed magnetic field intensity, Phi and Phi' are 
25 known and computed electrical potentials, lambda is a scaling factor to provide a 
weighting between the two objectives, and f is a general function which provides a 
measure of the difference between the two fields, for example absolute magnitude 
difference, difference in pattern, difference in relative timings. 

30 Distributed source inverse algorithms are commonly used in solving the inverse problem 
of electrocardiography, although point source inverses are also used to a lesser degree. 
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There are two main source formulations, that of a potential-based inverse and that of 
activation time inverses. The potential-based formulation determines a temporally 
varying electrical potential field for a given surface, while an activation time inverse 
algorithm determines the time at which the electrical potential wave front passes each 
5 point in space. 

Having obtained magnetic recordings from a subject and devised an accurate geometric 
model of the subject, the invention permits localisation of the source of the magnetic 
field. Magnetic fields are typically recorded at between 20 and 100 sites and provide a 
10 combination of gradient magnetic vector field and absolute magnetic field vector 
recordings. 

Figure 10 illustrates magnetic field vectors 1000 that are recorded from a human subject 
shown at 1002 and 1004 using a SQUID 25 seconds apart. The magnetic fields represent 
15 a summation of the electrical activity occurring within the stomach indicated generally at 
1006. Images 1002 and 1004 show 37 vector magnetic/gradient field recordings at 20 
physical locations just above the body surface. One such physical location is indicated at 
1008. At five of the locations, there are multiple channels recording at different vector 
directions, for example the location indicated at 1010. 

20 

Using simulation studies and under controlled conditions, it is possible to localise the 
source of electrical activity to less than 1mm. The localisation is performed using both 
electrical recordings and magnetic recordings, and a combination of the two as the 
"measured" field. 

25 

Figure 11 illustrates point source localisation 1100 of a site of focal activity in the 
stomach below the fundus. Recording electrodes are shown in 1102, for example at 1104 
and magnetic sensors near the skin surface indicated at 1106. 
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In diagram 1108, magnetic fields are used to determine the site of focal activity on the 
stomach. In this case the dipole centre has been localised to within 1mm under controlled 
conditions. 

5 The model defined in accordance with the invention has the potential to provide valuable 
insight to activity which cannot be easily observed or measured using existing imaging 
methods. Validation of the model can be conducted using analytic models or validation 
experiments. Analytic models usually involve simplified geometries and known initial 
and boundary conditions are applied to obtain a known solution. A more thorough test 
10 involves the use of validation experiments where measurements are obtained both 
internal and external of the subject. In this way, the computed solutions can be compared 
to those measured directly within the torso. 

One of the areas of application for the invention is in the modelling of gastroparesis and 
15 gastric uncoupling. One type of gastroparesis can occur when nerves to the stomach are 
damaged or stop working. The vagus nerve controls the movement of food through the 
digestive tract. If the vagus nerve is damaged, the muscles of the stomach and intestines 
do not work normally, and the movement of food is slowed or stopped. 

20 Gastroparesis is relatively prevalent in patients with type 1 diabetes and is detrimental to 
the patient's health as the delay in gastric emptying inhibits the control of blood glucose 
levels. At least 20% of people with type 1 diabetes develop gastroparesis. This condition 
also occurs in people with type 2 diabetes although with less frequency. 

25 The cause of gastric uncoupling can be mechanical where external intervention, for 
example surgery or trauma, creates an electrical break in the stomach wall. The normal 
pacemaking wave is no longer able to propagate past the site of the break and pacemaker 
cells distal to the break begin pacing the distal stomach at a slower rate than the normal 
dominant frequency. 

30 
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Using the invention, a model is created of a slice through the stomach wall along the 
greater curvature. The model includes Interstitial Cells of Cajal (ICCs) along with 
circular and longitudinal smooth muscle layers. 

5 The model is first executed in the absence of any mechanical uncoupling. In this 
situation, the distal regions of the tissue model were entrained at the frequency of the 
dominant pacemaker which was approximately 3.06 cpm (0.051Hz). Transmembrane 
potentials and the power spectrum of the electrical control activity (EC A) over time are 
displayed in the signal output window in Figure 12. 

10 

The sample output 1200 shows at 1202 the computed transmembrane potentials in mV 
and time in seconds and shows at 1204 the power spectrum of the top trace where the 
frequency is in Hz. 

15 A conduction blockage was then introduced 50% of the way from the dominant 
pacemaker site to the terminal antrum. Proximal to the blockage, the dominant ECA 
frequency remained unchanged, but distal to the blockage the dominant frequency was 
reduced to 2.22 cpm (0.038Hz) as shown in Figure 13. 

20 Figure 13 shows at 1300 sample output from the distal region of the tissue model after 
mechanical uncoupling. Graph 1302 shows the transmembrane potential in mV plotted 
against time in seconds and 1304 shows the power spectrum of the signal trace with a 
decreased dominant frequency shown in Hz. 

25 This bradygastria due to mechanical uncoupling coincides well with measured 
bradygastrias of the same origin. 

The examples described above with reference to Figure 3 to 12 have particular 
application in data interpretation relating to the gastrointestinal tract. It will be 
30 appreciated that the same techniques could be applied to data interpretation related to 
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gastric, intestine, and cardiac activity, and could be used to model and analyse any organ, 
combination of organs, or part(s) thereof. 

The foregoing describes the invention including preferred forms thereof. Alterations and 
5 modifications as will be obvious to those skilled in the art are intended to be incorporated 
within the scope hereof, as defined by the accompanying claims. 
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